Mirror quiescence and high-sensitivity position measurements with feedback 
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We present a detailed study of how phase-sensitive feedback schemes can be used to improve 
the performance of optomechanical devices. Considering the case of a cavity mode coupled to 
an oscillating mirror by the radiation pressure, we show how feedback can be used to reduce the 
position noise spectrum of the mirror, cool it to its quantum ground state, or achieve position 
squeezing. Then, we show that even though feedback is not able to improve the sensitivity of 
stationary position spectral measurements, it is possible to design a nonstationary strategy able to 
I increase this sensitivity. 
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I. INTRODUCTION 



Mirrors play a crucial role in a variety of precision measurements like gravitational wave detection and atomic 
force microscopes [Q. In these applications one needs a very high resolution for position measurements and a good 
K*" " control of the various noise sources, because one has to detect the effect of a very weak force [||J^]. As shown by 
the pioneering work of Braginsky even though all classical noise sources had been minimized, the detection of 
gravitational waves would be ultimately determined by quantum fluctuations and the Heisenberg uncertainty principle. 
Quantum noise in interferometers has two fundamental sources, the photon shot noise of the laser beam, prevailing 
, at low laser intensity, and the fluctuations of the mirror position due to radiation pressure, which is proportional 
T-H ' to the incident laser power. This radiation pressure noise is the so-called "back-action noise" arising from the fact 
] that intensity fluctuations affect the momentum fluctuations of the mirror, which are then fed back into the position 

■ by the dynamics of the mirror. The two quantum noises are minimized at an optimal, intermediate, laser power, 
Q-i' yielding the so-called standard quantum limit (SQL) ^j^. Real devices constructed up to now are still far from the 
■fjj standard quantum limit because quantum noise is much smaller than that of classical origin, which is essentially given 
^ by thermal noise. In fact, present interferometric gravitational wave detectors are limited by the Brownian motion 
^ ■ of the suspended mirrors 0, which can be decomposed into suspension and internal (i.e. of internal acoustic modes) 

thermal noise. Therefore it is very important to establish the experimental limitations determined by the thermal 

■ noise, and recent experiments §,01 go in this direction. 

. !^ I Recently it has been reported the first experimental evidence of the reduction of thermal noise by means of the 

■ radiation pressure of an appropriately modulated laser light incident on the back of the mirror . The method was 
^ [ based on a phase-sensitive feedback control proposed in Ref. JT2| : detect the mirror displacement through a homodyne 
- - ' measurement, and then use the output photocurrent to realize a real-time reduction of the mirror fluctuations. The 

proposed scheme is a sort of continuous version of the stochastic cooling technique used in accelerators , because 
the feedback continuously "kicks" the mirror in order to put it in its equilibrium position. This proposal has been 
experimentally realized in Ref. using the "cold damping" technique JT4| , which amounts to applying a viscous 
feedback force to the oscillating mirror. In the experimental studies of optomechanical systems performed up to 
now, the effects of quantum noise are blurred by thermal noise and the experimental results can be well explained in 
classical terms (see for example |l5j). However, developing a fully quantum description of the system in the presence 
of feedback is of fundamental importance, for two main reasons. First of all it allows to establish the conditions under 
which the effects of quantum noise in optomechanical systems become visible and experimentally detectable. We 
have recently shown in Ref. [ p^ that there is an appreciable difference between the classical and quantum description 
of feedback already at liquid He temperatures. Moreover, a completely quantum treatment allows to establish the 
ultimate limits of the proposed feedback schemes, as for example, the possibility to reach ground state cooling of 
a mechanical, macroscopic degree of freedom. In Ref. [Q, a quantum treatment of stochastic cooling feedback has 
been already presented, based however on a master equation description which is not valid at very low temperatures 
A consistent quantum description of both stochastic cooling and cold damping feedback schemes, valid at all 
temperatures, has been presented in |ll6|] , and recently a discussion of the quantum limits of cold damping has been 
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presented in ||T^ . The present paper will extend and generalize the results of jT^,|l^ , allowing us to make a detailed 
comparison of the two feedback schemes, and to establish all their potential applications. In particular we shall see 
that both schemes can achieve ground state cooling of an oscillating mode of the mirror, and that, in an appropriate 
limit, the "stochastic cooling" feedback of Ref. [Q, can even break the standard quantum limit, achieving steady 
state position squeezing. The experimental realization of these quantum limits in optomechanical systems is extremely 
difficult, but the feedback methods described in this paper may be useful also for microelectromechanical systems, 
where the search for quantum effects in mechanical systems is also very active . 

Thermal noise reduction is important, but is not the only relevant aspect. What is more important, expecially for 
gravitational wave detection 0|, or for metrology applications |^0|, is to improve the sensitivity, i.e., the signal to 
noise ratio (SNR) of position measurements Q. Both the stochastic cooling scheme of Ref. and the cold damping 
scheme of Ref. pT| cool the mirror by overdamping it, thereby strongly decreasing its mechanical susceptibility at 
resonance. Cooling is therefore achieved through the suppression of the resonance peak in the noise power spectrum. 
This suggests that both feedback schemes cannot be directly applied to improve the sensitivity for the detection of 
weak forces, because the strong reduction of the mechanical susceptibility at resonance means that the mirror does 
not respond both to the noise and to the signal. We shall see that this is true only in stationary conditions, i.e., we 
shall prove that the stationary spectral SNR is never improved by feedback. However, as we have recently shown in 
, it is possible to use feedback with an appropriate nonstationary strategy, able to increase significantly the SNR 
for the detection of impulsive classical forces acting on the oscillator. Here we shall extend the results of by 
adopting a general description of nonstationary spectral measurements. 

The outline of the paper is as follows. In Sec. II we describe the model and derive the appropriate quantum 
Langevin equations. In Sec. HI we describe the stochastic cooling feedback scheme of Ref. [|l2| and the cold damping 
feedback using the quantum Langevin theory developed in pl| , p2| , and we make a detailed comparison of the two 
schemes. In Section IV we analyze the stationary state of the oscillating mirror, and we determine the conditions 
under which feedback can be used to achieve ground state cooling or position squeezing. In Section V we present a 
general description of nonstationary spectral measurement and we discuss the stationary limit in particular. Section 
VI describes how the sensitivity of position measurements can be improved by using feedback in a nonstationary way, 
and Section VII is for concluding remarks. 



II. THE MODEL 



The system studied in the present paper consists of a coherently driven optical cavity with a moving mirror 
(Fig. |l|) . This opto-mechanical system can represent one arm of an interferometer able to detect weak forces as those 
associated with gravitational waves jl] , or an atomic force microscope [|| . The detection of very weak forces requires 
having quantum limited devices, whose sensitivity is ultimately determined by the quantum fluctuations. For this 
reason we shall describe the mirror as a single quantum mechanical harmonic oscillator with mass m and frequency 
(jOra- Experimentally, the mirror motion is the result of the excitation of many vibrational modes, including internal 
acoustic modes. The description of the mirror as a single oscillator is however a good approximation when frequencies 
are limited to a bandwidth including a single mechanical resonance, by using for example a bandpass filter in the 
detection loop p3|. 
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FIG. 1. Schematic description of the system. The cavity mode is driven by the laser which, thanks to the beam splitter, 
provides also the local oscillator for the homodyne measurement. The signal is then fed back to the mirror motion. 



The optomechanical coupling between the mirror and the cavity field is realized by the radiation pressure. The 
electromagnetic field exerts a force on the movable mirror which is proportional to the intensity of the field, which, at 
the same time, is phase-shifted by 2kq, where k is the wave vector and q is the mirror displacement from the equilibrium 
position. In the adiabatic limit in which the mirror frequency is much smaller than the cavity free spectral range 
c/2L {L is the cavity length) ||2j], one can focus on one cavity mode only because photon scattering into other modes 
can be neglected, and one has the following Hamiltonian pS] 



H = fiujjj'^b + nuj,n {P^ + Q^) - TiGb^bQ + ifiE {b'^ e 



be 



(1) 



where b is the cavity mode annihilation operator with optical frequency uJc^ and E describes the coherent input field 
with frequency ~ driving the cavity. Moreover, Q and P are the dimcnsionless position and momentum operator 
of the movable mirror, with [Qj-P] = V^j and G — {wc/ L)yJh/2mijJm is the coupling constant. Since we shall focus 
on the quantum and thermal noise of the system, wc shall neglect all the technical sources of noise, i.e., we shall 
assume that the driving laser is stabilized in intensity and frequency. This means neglecting all the fluctuations of the 
complex parameter E. Including these supplementary noise sources is however quite straightforward and a detailed 
calculation of their effect is shown in Ref. ||2^ . Moreover recent experiments have shown that classical laser noise 
can be made negligible in the relevant frequency range The adiabatic regime ujm ^ c/2L we have assumed in 

Eq. implies w„i ^ and therefore the generation of photons due to the Casimir effect, and also retardation and 
Doppler effects are completely negligible. 

The dynamics of the system is not only determined by the Hamiltonian interaction (|l|), but also by the dissipative 
interaction with external degrees of freedom. The cavity mode is damped due to the photon leakage through the 
mirrors which couple the cavity mode with the continuum of the outside electromagnetic modes. For simplicity we 
assume that the movable mirror has perfect reflectivity and that transmission takes place through the other, "fixed" , 
mirror only. We indicate the photon decay rate at the fixed mirror by 7c. Then, the quantity E is related to the 
input laser power p by _E = p7c/?iwo- The mechanical oscillator, which may represent not only the center-of-mass 
degree of freedom of the mirror, but also a torsional degree of freedom as in or an internal acoustic mode as in 
undergoes Brownian motion caused by the uncontrolled coupling with other internal and external modes at thermal 
equilibrium. 

The dynamics of the system can be described by the following set of coupled quantum Langevin equations (QLE) 
(in the interaction picture with respect to TiLo^b'^b) 



Q{t) 
P{t) 

bit) 



^ i^rnP{t) , 

: -LJ„,Q{t) + W{t) - J^P{t) + Gb\t)b{t) , 

: - {iuc - iujQ + y ) b{t) + 2iGQ{t)b{t) + E 



(2a) 
(2b) 

(2c) 



where bin{t) is the input noise operator p7[ associated with the vacuum fluctuations of the continuum of modes 
outside the cavity, having the following correlation functions 
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{hn{t)bUt')) = {bl{t)b„,{t')) = , (3a) 

{hnit)blit'))=Sit^t'). (3b) 

Furthermore, VV(t) is the Brownian noise operator defined consistently with quantum mechanics p^ . It has the 
foUowing correlation functions 

mmit')) = ^^{Mt - 1') + t')} , (4) 

where 

J^rit) — / dto uj cos{ujt) coth I — — — I , (5a) 







2kBT^ 

^i{t) = — I duj a;sin(a;i) , (5b) 





with T the bath temperature, "fm the mechanical decay rate, ks the Boltzmann constant, and w the frequency 
cutoff of the reservoir spectrum. The antisymmetric part, J-i, of Eq. (^, is a direct consequence of the commutation 
relations for the Brownian noise operator, and the symmetric part, explicitely depends on temperature and becomes 
proportional to a Dirac delta function when the high temperature limit fc^T 3> hvo first, and the infinite frequency 
cutoff limit — !■ oo later, are taken. Eqs. (^ and (^) show the non-Markovian nature of quantum Brownian motion, 
which becomes particularly evident in the low temperature limit | |2^ , p9| . Therefore, the exact QLE (|^) reduce to the 
standard ones p^ ] in the limit w oo. It is also important to stress that the quantum Langevin description of 
quantum Brownian motion given by Eq. (|^) is more general than that associated with a master equation approach, 
because it is valid at all temperatures and it docs not need any high temperature limit [0. 

In standard interferometric applications, the driving field is very intense. Under this condition, the system is char- 
acterized by a semiclassical steady state with the internal cavity mode in a coherent state |/3), and a new equilibrium 
position for the mirror, displaced by G\P\'^ /uj„i with respect to that with no driving field. The steady state amplitude 
is given by the solution of the nonlinear equation 

P = , (6) 

f +^^c-i^o-2z£|/3|2' 

which is obtained by taking the expectation values of Eqs. (||), factorizing them and setting all the time derivatives to 
zero. Eq. (^) shows a bistable behaviour which has been experimentally observed in [ pO[ . Under these semiclassical 
conditions, the dynamics is well described by linearizing the QLE (|^) around the steady state. If we now rename with 
Q(t) and b{t) the operators describing the quantum fluctuations around the classical steady state, we get 

Q{t) = uj^P{t) , (7a) 
P{t) - -c^„Q(t) - -i„,P{t) + G/3 [b[t) + b\t)\ + W{t) , (7b) 

m = - (y + bit) + 2iGf3Q{t) + ^An{t) , (7c) 
where we have chosen the phase of the cavity mode field so that /3 is real and 

A = - tjo - — /3^ (8) 

is the cavity mode detuning. We shall consider from now on A = 0, which corresponds to the most common 
experimental situation, and which can always be achieved by appropriately adjusting the driving field frequency luq. 
In this case the dynamics becomes simpler, and, introducing the field phase quadrature Y{t) = i (6^(i) — 6(t)) /2 and 
field amplitude quadrature X{t) — (6(i) + b^t)) /2, one has that only the phase quadrature Y{t) is affected by the 
mirror position fluctuations Q{t), while the amplitude field quadrature X{t) is not. In fact, the hnearized QLE (|^) 
can be rewritten as 

Q{t) = uj,nP{t) , (9a) 
P{t) = -Lj„,Q{t) - 7„P(t) + 2Gl3X{t) + W{t) , (9b) 

Y{t) - -^Y{t) + 2G(3Q{t) + ^Yn{t) , (9c) 
i:(i) = -|x(i) + ^X„(t), (9d) 
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where we have introduced the phase input noise Yin(t) = i ^6j„(t) — bin{t)j and the amphtude input noise Xin(t) — 

III. POSITION MEASUREMENT AND FEEDBACK 

Usuahy the movable mirror is used as a ponderomotive meter to detect smah forces acting on it Q. Thus, we 
introduce an additional Hamiltonian term describing the action of a classical external force f{t), that is 

Hext = -Qfit) . (10) 

Information about such a force can be obtained by looking at the mechanical oscillator position Q(t). The position 
measurement is commonly performed in the large cavity bandwidth limit 7c 3> G/3, Wm, when the cavity mode 
dynamics adiabatically follows that of the movable mirror and it can be eliminated, that is, from Eq. 



(11, 



and X{t) ~ Xin{t) / ^/j^ from Eq. (pq). Performing a continuous homodyne measurement of the phase quadrature 
Y{t) means therefore continuously monitoring the real time dynamics of the oscillator position Q(t), which, in turn, 
implies detecting the effects of classical force f(t). The experimentally detected quantity is the output homodyne 
photocurrent [|l]j22|j3|] 



Youtit) = 2ij^eY{t) - VvyZit) , (12) 

where 77 is the detection efficiency and Yll^{t) is a generalized phase input noise, coinciding with the input noise Yin(t) 
in the case of perfect detection rj — 1, and taking into account the additional noise due to the inefficient detection 
in the general case 77 < 1 p2| . This generalized phase input noise can be written in terms of a generalized input 
noise br,{t) as YZ{t) = i [6jj(t) — bri{t)\. The quantum noise 6,,(t) is correlated with the input noise 6m (t) and it is 
characterized by the following correlation functions 

{b,{t)b,{t')) = {bl{t)b,{t'))^0, (13a) 
{b^it)bl{t'))=d{t-t'), (13b) 
{b^nit)bl{t')) = {b,it)bl{t')) ^ ^j6it - t'). (13c) 

The output of the homodyne measurement may be used to devise a phase-sensitive feedback loop to control the 
dynamics of the mirror. For example, we have proposed in Ref. [jl^ to reduce the effects of thermal noise on the 
mirror by feeding back the output homodyne photocurrent in an appropriate way. The proposed scheme is a sort of 
continuous version of the stochastic cooling technique used in accelerators , because the homodyne measurement 
provides a continuous monitoring of the oscillator's position, and the feedback continuously "kicks" the mirror in order 



to put it in its ec 
realized in Ref. 



uilibrium position. Our proposal of cooling the mirror using a feedback loop has been experimentally 
rT] (see also jl^), using a different method, the so-called "cold damping" technique [Q. This latter 
feedback scheme shares some analogies with that proposed in Ref. [|2| and amounts to applying a viscous feedback 
force to the oscillating mirror. In the experiment of Refs. |]lT| , p^ , the viscous force is provided by the radiation 
pressure of another laser beam, intensity-modulated by the time derivative of the homodyne signal. 

The effect of the feedback loop has been described using quantum trajectory theory and the master equation 
formalism in Ref. and a classical description neglecting all quantum fluctuations in Ref. |l^Jl^. Here we shall 
use a more general description of feedback based on QLEs for Heisenberg operators, first developed in Ref. |^ and 
generalized to the non-ideal detection case in Ref. (see also ||3^ for a comparison between these quantum feedback 
approaches and general quantum control theories). This general quantum description of feedback will allow us to 
compare the two different feedback schemes, the stochastic cooling scheme of Ref. [|l2|, and the cold damping scheme 
of Ref. ||ll|,|l^. Moreover, the present quantum treatment will allow us to show that in the presence of feedback 
the radiation quantum noise has important effects, and that a classical stochastic treatment of the dynamics of the 
system is generally inadequate. Our treatment explicitely includes the limitations due to the quantum efficiency of 
the detection, but neglects other possible technical imperfections of the feedback loop, as for example the electronic 
noise of the feedback loop (discussed in ) , or the fluctuations of the laser beam used for the feedback in the cold 
damping scheme. 



5 



A. Stochastic cooling 



Let us first consider the stochastic cooling scheme of Ref . ||T^ . In this scheme, the feedback loop induces a continuous 
position shift controlled by the output homodyne photocurrent Yout (t) ■ This effect of feedback manifests itself in an 
additional term in the QLE for a generic operator 0{t) given by [EJ 



Ofb{t) = i^Yout{t~T) [g,,P{t),0{t)] , 



(14) 



where t is the feedback loop delay time, and gsc is a dimensionless feedback gain factor. The feedback delay-time is 
essentially determined by the electronics involved in the feedback loop and is always much smaller than the typical 
timescale of the mirror dynamics. It is therefore common to consider the zero delay-time limit, r — > 0. This limit is 
however quite delicate in general j2^j2^. In fact, Youtit — T), being an output operator, commutes with [gscP{t): 0{t)] 
for any nonzero r, but this is no more true when t — Q. Therefore, one has to be careful with ordering in the zero 
delay-time limit. However, with the choice of Eq. ( p^ for the feedback term, the only nonzero commutator in the 
QLE of Eq. (^) is [gsc^'(Oj which, being a c-number, does not create any ordering ambiguity. Therefore one 

has exactly the same equations one would have by putting directly r = in Eq. that is. 



Q{t) = oj„,P{t) + 9sclcY{t) - ^-f^l^Ylit) , 

P{t) = -Lu^Qit) - -,^P{t) + 2GpX{t) + W{t) + fit) 

Y{t) = -^Y{t) + 2G(3Q{t) + ^y.n{t) , 



X(t) 



X{t) 



(15a) 
(15b) 
(15c) 

(15d) 



where we have used Eq. (^^. After the adiabatic elimination of the radiation mode (see Eq. (^_1|)), the above equations 
reduce to 



Q{t) = uj^P{t) + AGfigscQit) + ^c9scYUt) - -^J-Klit) , 

z y 77 

Pit) = -u^Qit) - 7,„P(t) + ^X„,(i) + Wit) + fit) . 



lie 



(16a) 
(16b) 



The solution of these QLE for the conjugate operators Q(t) and Pit) can be easily obtained by performing the Laplace 
transform, and they will be useful in the following. Their expression is 



Q(t) = ifQ(i)Q(0) + Xsc(i)P(0)+ / dt'xscit^)fit-t!) 



dt'Kgit') 



fTc9s.YUt^t')-^-f^Ylit-t') 



+ / dt'xscit') 



wit-t') + ^x,nit-t') 



Pit) = Kpit)PiQ) - Xscit)QiQ) + / dt' Kpit')fit - t') 



dt' Kpit') 



Wit^t') + ^X^nit-t') 



dt' Xscit') VTcgscYn it-t')-^ ^Yg it - t') 



(17a) 



(17b) 



We have introduced the time-dependent susceptibility Xscit) describing the response of the movable mirror in the 
presence of the stochastic cooling feedback 



Xscit) 



\/^m -7m (^) 



^e-(i+9i)7™*/2sin 



1-51 



(18) 
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and the two related response functions 

Kgit) = (19) 

We have also rescaled the feedback gain and defined gi — —'iGfigschm- 

B. Cold damping 

Cold damping techniques, that is, the possibility to use a feedback loop to reduce the effective temperature of a 
system well below the operating temperature, have been applied in classical electromechanical systems for many years 
and only recently they have been proposed to improve cooling and sensitivity at the quantum level js^. This 
technique is based on the application of a negative derivative feedback, which increases the damping of the system 
without correspondingly increasing the thermal noise QjsJ] . This technique has been succcsfuUy applied for the first 
time to an optomechanical system composed of a high-finesse cavity with a movable mirror in the experiments of 
Refs. In these experiments, the displacement of the mirror is measured with very high sensitivity ||], and the 

obtained information is fed back to the mirror via the radiation pressure of another, intensity-modulated, laser beam, 
incident on the back of the mirror. Cold damping is obtained by modulating with the time derivative of the homodyne 
signal, in such a way that the radiation pressure force is proportional to the mirror velocity. The servo-control force 
then corresponds to a viscous force. The results of Refs. pl] , p^ referred to a room temperature experiment, and have 
been explained using a classical description. The quantum description of cold damping in this optomechanical system 
has been presented in |l^] (see also Ref. |^^)), and we shall follow this treatment. 

In the quantum Langevin description, cold damping feedback scheme implies the following additional term in the 
QLE for a generic operator 0{t) p^ , 

Ofbit) = ^—Y,^t{t^T) [g,dQ{t),0{t)] . (21) 

As for the stochastic cooling feedback case, one has only one nonzero feedback term in the QLE of the system (||), 
which in this case is [gcdQ{t), Pit)]- Since this commutator is a c-number also in this case, we do not have any ordering 
problem in the zero delay-time limit, and the QLE for the cold damping feedback scheme become 

Q{t) = uj,nPit) , (22a) 
Pit) = -LUraQit) - jmPit) + 2G/3X(t) - g.dYit) + -^Y^t) + Wit) + fit) , (22b) 



Yit) - -^Yit) + 2GPQit) + ^r.„(t) , (22c) 

7c , 



X{t) = -fX{t) + ^XUt)- (22d) 
Adiabatically eliminating the cavity mode, one has 

Qit)=Lu„,Pit), (23a) 

Pit) = -COrnQit) ~ InrPit) + ^X„(0 + Wit) + fit) - ^^^Q(0 - -^^^(i) + -^F^Jt). (23b) 



Notice that the modulation with the derivative of the homodyne photocurrent implies the introduction of two new 
quantum input noises, Yinit) and Y^it), whose correlation functions can be simply obtained by differentiating the 
corresponding correlation functions of Yin it) and Y^j-^ it) . We have therefore 

{Ynit)Ynit')) ^ {Ynit')Ynit)) = {Ylit)Ylit')) = {Ylit')Ylit)) ^ -^t - t!), (24a) 
iYlit)YM = (YUmiit)) = -Vv'Sit t'). (24b) 
{X.nit)Y:iit')) = -{Ylit')XUt)) = -^^^5it - t'). (24c) 
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In this case the solution of the adiabatic QLE reads 



Q{t) = Kit)QiO) + Xcd{t)P{Q) + / dt' xcd{t')f{t - t') 

Jo 

+ I dt'xcd{t~t') 



(25) 



and P{t) = Q{t)/u!jrn where we have introduced the time-dependent susceptibiUty in the case of the cold damping 
feedback scheme 



Xcd{t) 



^g-(l+92)7™t/2gijj 



,,2 f 1+92 ^ 



'~fm V 2 

and the related response function 



(26) 



K{t) = l-w™ / dt'xcdit') 



We have again rescaled the feedback gain and defined g2 ~ '^Gj3ujmgcdllmlc- 



C. Comparison between the two feedback schemes 



(27) 



The two sets of QLE for the mirror Heisenberg operators, Eqs. and (|23D, show that the two feedback schemes 
are not exactly equivalent. They are however physically analogous, as it can be seen, for example, by looking at the 
differential equation for the displacement operator Q{t). In fact, from Eqs. (|l^) one gets 



Q{t) + (1 + gi) j.mQ{t) + {ujI + -flgi) Qit) = 



2G/3 



X,„it)+Wit)+f{t) 



y^gscYinit) 



for the stochastic cooling scheme, while from Eqs. (E3[) one gets 



Q{t) + (1 + 52) 7mQ(t) + ujl,Q{t) = w„ 



2G/3 



/7c 



X^n{t)+Wit)+f{t) 



(28) 



(29) 



for the cold damping scheme. These equations shows that in both schemes the main effect of feedback is the modifi- 
cation of mechanical damping 7^ 7m (1 + gi) (* = 1, 2). In the stochastic cooling scheme one has also a frequency 
renormalization -I- 7,^31 , which is however usually negligible since the mechanical quality factor Q = i^mj^m 

is always large. Moreover, in the two cases the position dynamics is affected by similar, even though not identical, 
noise terms. This comparison shows that the stochastic cooling scheme of Ref. is also able to provide a cold 
damping effect of increased damping without an increased temperature [16[. 



IV. STATIONARY STATE AND COOLING 



We now study the stationary state of the movable mirror in the presence of the two feedback schemes, which is 
obtained by considering the dynamics in the asymptotic limit i — > cx3. We shall see that both feedback schemes 
are able to lower the effective temperature of the system, and that, in particular limits, the steady state can have 
interesting quantum features. In fact, both schemes are able to achieve ground state cooling, and the stochastic 
cooling feedback is even able to achieve steady state position squeezing. 
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A. Stochastic cooling feedback 



Using the solution (17a), one has 



dt' dt"KQ{t')KQ{t")c,{t' ~t")+ dt' dt"xsc{t')Xsc{t")c2{t' -t"), (30) 



{Q^)st = lim {Q{tf) 



where ci{t) is the stationary symmetrized correlation function of the noise term ni{t) = ^/j^gscYmit) — ^^/^^/^(i)j 
C2(t) is the stationary symmetrized correlation function of the noise term n2{t) ~ VV(t) 



2G/3 



Xin{t), and we have 



used the fact that rii (t) and n2 {t) are uncorrelated. Using the correlation functions (||) , (||) , and (|T3|) , one gets 



ciit) 



64r7G2/32 



m 



C2{t) = '—5{t) 



-Tr{t). 



7c 27rw,i 

The expression for {Q'^)st is obtained using Eqs. (|T^), (|l9|), and ( |3l| ) in Eq. 



7c7m5? 



1 + 



.91 



2fl2 



1287?G2^2 (l+5i)(Q'+.9i) 7c7m (l+.9i)(Q'+5i) 



{Q'^)bm- 



(31a) 
(31b) 

(32) 



The term {Q'^)bm is the contribution of the mirror quantum Brownian motion, whose general expression is obtained 



by rewriting J-r{t' — t") in Eq. (|30|) in terms of its Fourier transform J-'riw) (see Eq. 



to get 



where 



{Q'^)bm 



dw 



In 



2tt 2uJr, 



) coth 



/ huj 



\2kBT 



lx^c(w)l^ 



^li + 5i7m - + «W7™ (1 + gi) 



(33) 



(34) 



is the frequency-dependent susceptibility of the mirror in the stochastic cooling feedback scheme. The general analyti- 
cal expression of the quantum Brownian motion term (Q^) em, valid in any range of parameters, is cumbersome and has 
been obtained in ||2^,^. However, in typical optomechanical experiments [^,^,^,^ it is always fi^m ^ h^^m ^ ksT, 
and it is possible to see [2^ ] that, in this limiting case, the classical approximation coth(?ia;/2fcBT) ~ 2kBT /huj (which 
is equivalent to approximate J-r{t) ~ {'^rnksT /Tiujm) S{t)) can be safely used in (|33|), so to get 



{Q'^)bm 



knT 



2hu;,n (l + 5i)(Q^+5i) 



Finally it is 



{Q^)st = 



9l 



1 



.91 



2huJr, 



(l+5i)(Q'+5i) 



8ryC(i + .9i)(Q' + .9i) 

where we have introduced the rescaled, dimensionless, input power of the driving laser 

16G2/32 64G2 



c 



(35) 



(36) 



(37) 



Eq. (^) coincides with the corresponding one obtained in ||T^ using a Master equation description of the stochastic 

cooling feedback scheme. 

An analogous procedure can be followed to get the stationary value (P^)st. Using Eqs. ( |l7b|) , (^8|), (20), (|3l|), and 
(^), one obtains the general expression 



.9? 



■9i 



8??C(i + 5i)(Q2 + 5i) s a + gi) iQ^ + gi) 



BM, 



(38) 
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where the quantum Brownian motion contribution is now given by 



{P')bm 



duo 7ni 

27r 2w, 



-w coth 



2kBT 



\Xsc{(^)\' 



(39) 



In this case, the classical, high-temperature, approximation colh{fiU) /2kBT) ~ 2kBT/l%u> has to be made with care, 
because, due to the presence of the term, the integral ( |39| ) has an ultraviolet divergence in the usually considered 
zu ^ oo limit (see also Eq. (jsj)). This means that, differently from {Q'^)bMi the classical approximation for {P'^)bm 
is valid only under the stronger condition fivo <^ kBT and that in the intermediate temperature range hvj 3> 
fcsT 3> fiujm (which may be of interest for optomechanical systems), one has a correction of order In {h,w /ksT). One 
has therefore |2al 



(P2 



kuT 



BM — 



9l 
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2hiOrn (l+5l)(S'+5l) 



■In 



\2TTkTj 



so that one finally gets 



(P')st = 



9l 



8'7C(l + 5i)(Q'+.9i) 



C , kBT 



2flL0„ 



gl + Q^ + 91 

(l + 5i)(Q^+5i) 



■In 



\2TTkT 



(40) 



(41) 



This expression coincides with the corresponding one obtained in using a Master equation description, except for 
the logarithmic correction, which however, in the case of mirror with a good quality factor Q, is quite small, even in 
the intermediate temperature range hzu ^ kBT ^ fu-Om- 

A peculiar aspect of the stochastic cooling feedback scheme, which has not been underlined in jl^], is its capability 
of inducing steady-state correlations between the position and the momentum of the mirror, i.e., the fact t hat { QP + 
PQ)st 7^ 0. This correlation can be evaluated in the same way as above, starting from Eqs. (17a) and ( [l7b| ), and 
getting 



{QP + PQ)s 



dt' dt"KQ{t')Xsc{t")c^{t' -t")+ dt' dt"Kp{t')Xsc{t")c2(t' ~t") 



(42) 



Then, using Eqs. (^8|), ([l9|), (po|), and (|3l|), and performing the classical approximation on the quantum Brownian 
motion contribution (there is no ultraviolet divergence for — s- oo in this case), one gets 



{QP + PQ)s 



knT 



giQ 



9l 



Q 



C 

8 ' 2hL0ra J {I + gi)iQ'' + gi) 8?7C(l + <?l)(Q2 + gi) 



(43) 



Each steady state expression (36), ( ^l]) and ( ^3|) has three contributions: the thermal term due to the mirror Brownian 
motion, the back action of the radiation pressure, proportional to the input power and the feedback-induced noise 
term proportional to and inversely proportional to the input power. At sufficiently large temperatures, the thermal 
noise contribution is much larger than the others and the mirror dynamics is faithfully described in terms of classical 
stochastic equations. This classical description amounts to neglect all the radiation input noises into the evolution 
equations of the system, so that W(t) is the only noise acting on the system. This classical description has been 
succesfully used in Rcfs. ||ll|,|l^ to account for the experimental data, in the case of a cold damping feedback scheme 
at room temperature. It is however evident that the radiation back action and the feedback-induced noise cannot 
be neglected in general. For example, the classical approximation for (Q^)^* suggests that it would be possible to 
localize the mirror without limit, i.e., {Q'^)st ^ 0, using an ever increasing feedback gain gi and keeping the input 
power fixed, while this is no more true as soon as the feedback-induced noise term proportional to g^ is included. 

The stochastic cooling feedback scheme has been introduced in | p2[ as a promising method for significantly cool- 
ing the cavity mirror. Let us therefore consider the optimal conditions for cooling, and the cooling limits of this 
scheme. The interesting quantity is the stationary oscillator energy Ust, which, neglecting the logarithmic correction 
of Eq. (|4l|)), can be written as 



{P')st] 



huJri 



gi 



11 + 2Q2+5,) 



<(l + 5i)(Q'+5i) 



+ C + 



^kBT\ (.g? + 2Q2 



flUOr, 
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(i + gi)(Q2 + gi) 



(44) 



It is evident from Eq. (^J) that the effective temperature is decreased only if both Q and gi are very large. At the 
same time, the additional terms due to the feedback-induced noise and the back-action noise have to remain bounded 
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for Q — > oo and gi — s- oo, and this can be obtained by minimizing Ust with respect to ^ keeping Q and gi fixed 
(physically this means optimizing the input power p at given gi and Q) . It is possible to check that these additional 
terms are bounded only for very large Q, that is, if Q/gi — > oo and in this case the minimizing rescaled input power 
is Copt — 9i/\/V- Under these conditions, the steady state oscillator energy becomes 



huJr, 



1 2kBT 1 



huJr, 
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(45) 



showing that, in the ideal limit rj = 1, gi oo, C ^ 9i ~^ oo, Q/gi oo, the stochastic cooling feedback scheme is 
able to reach the quantum limit Ust — f'-'-^m/'^, i-e., it is able to cool the mirror down to its quantum ground state. The 
behavior of the steady-state energy is shown in Figs. and H, where Ust (in zero-point energy units fnoml'^) is plotted 
as a function of the rescaled input power Q. In Fig. g, 2Ust/^i^m is plotted for increasing values of 171 (a: gi = 10, b: 
gi = 10^, c: gi = 10^ d: gi = lO'^) at fixed Q = lO'^, and with kBT/Tiuj^ = 10^ and i] = 0.8. The figure shows the 
corresponding increase of the optimal input power minimizing the energy, and that for high gain values, ground state 
cooling can be essentially achieved, even with a nonunit detection efficiency. In Fig. H, 2Ustlf>^m is instead plotted 
for increasing values of the mechanical quality factor Q (a: Q = 10'^, h: Q— 10^, c: Q— 10^) at fixed gi = 10''. The 
figure clearly shows the importance of Q in stochastic cooling feedback and that ground state cooling is achieved only 
when Q is sufficiently large. 




FIG. 2. Rescaled steady-state energy 2Ust/'hujm versus the rescaled input power C,, plotted for different values of gi (a: 
gi = 10, b: gi = 10^, c: gi = 10^, d: gi = lO'') at fixed Q = lO'^, and with ksT/naJm = 10^ and r; = 0.8. The optimal input 
power (^opt correspondingly increases, and for high gain values, ground state cooling can be achieved. 




FIG. 3. Rescaled steady-state energy 2Ust/TuVm versus 1^ for increasing values of the mechanical quality factor Q (a: Q — 10'', 
b: Q = lO'^, c: Q = lO'^) at fixed gi = lO'', and with kBT/hujm = 10^ and rj = 0.8. 

The possibility to reach ground state cooling of a macroscopic mirror using the feedback scheme of Ref. ||l2| was 
first pointed out, using an approximate treatment, in [3^ ], where the need of a very large mechanical quality factor is 
underlined. Here we confirm this result using the more general QLE approach. 
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The steady state of the mirror mode in the presence of stochastic cooling feedback shows other peculiar aspects and 
interesting limiting cases. Thanks to the linearization of the problem (see Eqs. this steady state is a Gaussian 
state, which however is never exactly a thermal state because it is always {Q^)st 7^ {P'^)st and {QP + PQ)st 7^ 0. Its 
phase space contours are therefore ellipses, rotated by an angle (j) = (1/2) arctan [{QP + PQ)st/ {{Q'^)st — (-F'^)st)] 
with respect to the Q axis. The steady state becomes approximately a thermal state only in the limit of very large 
Q (and » gi), as it can be seen from Eqs. (|3^), ( |4l] ) and (p3|). This thermal state approaches the quantum 
ground state of the oscillating mirror when also the feedback gain and the input power become very large. There 
are however other interesting limits in which the stochastic cooling feedback steady state shows nonclassical features. 
For example, the Gaussian steady state becomes a contractive state, which has been shown to be able to break the 
standard quantum limit in p6[ , when {QP + PQ)st becomes negative, and this can be achieved at sufficiently large 
feedback gain, that is, when gi > ri(^ {( + AkBT/hum) (see Eq. (^)). Finally, stochastic cooling feedback can be 
used even to achieve steady state position squeezing, that is, to beat the standard quantum limit (Q^)si < 1/4. The 
strategy is similar to that followed for cooling. First of all one has to minimize {Q'^)st with respect to the input power 
C at fixed gi and Q, obtaining 



giQ./1 + W 



91 



(46) 



This quantity can become arbitrarily small in the limit of very large feedback gain, and provided that gi 3> Q^. 
That is, difi'ercntly from cooling, position squeezing is achieved in the limit gi — > 00 (implying Q — > c»), and there 
is no condition on the mechanical quality factor. Under this limiting conditions, (Q^)st goes to zero as g^ and, 
at the same time, {P^)st diverges as gf^^, so that, in this limit, the steady state for the stochastic cooling feedback 
approaches the position eigenstate with Q = 0, that is, the mirror tends to be perfectly localized at its equilibrium 
position. The possibility to beat the standard quantum limit for the position uncertainty is shown in Fig. ^, where 
(Q^)st is plotted versus ( for two different values of the feedback gain, gi = 10^ (dotted hue), and gi — 10^ (full line), 
with Q — 10^, kBT/hoJm = 10^, and 77 = 0.8. For the higher value of the feedback gain, the standard quantum limit 
{Q^)st — 1/4 (dashed line) is beaten in a range of values of the input power 
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FIG. 4. Steady state position variance (Q^)at versus ^ for two values of the feedback gain, gi = 10^ (dotted line), and 
gi — 10^ (full line). The dashed line denotes the standard quantum limit {Q'^)st ~ 1/4, while the other parameters are: 
Q = 10*, kBT/hujm = 10^ and tj = 0.8. 



B. Cold damping feedback 



Now we characterize the stationary state of the mirror in the presence of cold damping. This stationary state has 
been already studied using classical arguments in ||ll|,Jl5| , while the discussion of the cooling limits of cold damping 
in the quantum case has been recently presented in ||l^. Here we shall generalize the results of ||l^ to the case of 
nonideal quantum efficiency 77 < 1, and we shall compare the cooling capabilities of the two feedback schemes. 

Using the solution (p5|) for the time evolution, one has 



{Q' 



= Km {Qitf) = / 



dt' / dt"x,,{t')xcd{t")cit' - t"), 



(47) 
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where c{t) is the stationary symmetrized correlation function of the noise term n{t) = ^^Xin{t')+W(t') — -^^Yin{t') + 



^^^±- Y_^(t') appearing in Eq. (^. Using the correlation functions (||), (^), ([l3|), and (p^), one gets 

4G2/32 



c{t) 



7c 



aid 



(48) 



Since in the cold damping case it is P(t) = Q{t)/ujm, it is straightforward to derive from Eq. ( |47| ) the expressions for 
{P'^)st and {PQ + QP)st, which are given by 



{PQ + QP)st = — lim ^(Qltf) = 0, 

UJm t ^oo at 



ip' 



1 



(49) 
(50) 



These stationary expressions can be rewritten in terms of the Fourier transforms of the noise correlation functions, in 
the same way as we have done for the Brownian motion term in the preceding subsection. Using Eqs. (^), (p6|), (|3^), 
and (E^), one has 



TT \Xcd[OJ)\ 
ZTT 



52 ^ , - CO 



coth 



huj 

2kBT 



duj UJ^ 2 



eAa;(w) + 



■ coth 



flLO 

2kBT 



where 



Xcrf(w) 



CJ^ - Cj2 + iuj-frn (1 + 32) 



(51) 

(52) 
(53) 



is the frequency-dependent susceptibility of the mirror in the cold damping feedback scheme, and Qi{lu) is a "gate" 
function, equal to one within the interval / and equal to zero outside. Notice that we have introduced not only the 
gate function 0j_^^^](a;) for the thermal noise term, but also the gate function 0Alj(^) for the feedback-induced noise 
term. In fact, it is easy to see that a frequency cutoff for the feedback is needed to avoid an ultraviolet divergence in 
the expression for {P'^)st- Moreover, from an experimental point of view, any feedback loop is active only within a 
finite bandwidth, which in this case is given by Aw. 

We first evaluate {Q'^)st- The contribution of the feedback-induced term generally depends upon the value of 
the feedback bandwidth Aw. There are two relevant experimental situations: a narrow bandwidth containing the 
mechanical resonance peak, that is, 7m (1 -t- 32) < Aw < ujm (configuration used in Ref. pl] , p^ ), or a wide bandwidth 
with a very large high frequency cutoff zuft, ^ Wm,7m(l +52)- However, since the factor |xcd(w)p in Eq. (|5^) is 
highly peaked around the resonance frequency Wm, {Q^)st is practically independent of the feedback loop bandwidth, 
as soon as 7m(l -I- 52) < Aw. In fact, either in the narrow bandwidth case, when the spectrum can be approximated 
by the constant term g'^/'^vC^ or in the case of a very large cutoff frequency, when the w^ dependence is kept, one gets 
the same result for the feedback-induced contribution, because 



duj LO 2 

7^ T \Xcd{uj)\ 



dcj 2 1 

TT \Xcd[uJ)\ = 7— — r. 

2n 2^„i[l + g2) 



(54) 



For the Brownian motion contribution we have the same situation described in the stochastic cooling case: the 
exact expression is cumbersome p8| , but in the commonly met condition hujrn ksT, the classical approximation 
coth(?iw/2fcBT) ~ 2kBT/Tujj can be made, and using Eq. (|5^ ) for both the thermal and the back-action contribution, 
one finally gets 



9l 



877C 8 2?iw„ 



1 



1 + 52 



(55) 



Notice that the corresponding expression for the stochastic cooling feedback ( |3q ) coincides with Eq. ( p5| ) in the limit 
Q> 1,51- 

Differently from {Q'^)st, {P'^)st depends upon the feedback loop bandwidth. In fact, in the large bandwidth case, 
the integrand in Eq. (j5^) tends to a constant at large frequencies, and in the limit of a very large cutoff frequency 
vcfb, the feedback-induced contribution becomes 
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(56) 



In the narrow bandwidth case instead, approximating the noise spectrum with the constant term .g|/477C, and using 
again Eq. (|5^ ) within Eq. (^), one gets a feedback-induced noise term contribution identical to that of {Q^)st of 
Eq. (p5|), which is independent of the feedback bandwidth. 

A potential ultraviolet divergence and a dependence upon the frequency cutoff w is present also in the quantum 
Brownian motion term. In fact, as we have seen in the preceding subsection, the classical expression for the thermal 
contribution to (P^)st, holds only in the limit of very large temperatures, kgT 3> hm, while, in the intermediate 
temperature regime huj„i <C kgT <^ hm, one has an additional logarithmic correction, so to get 



.2\BM 
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2flUJm 1 
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(57) 



Finally, the back-action term is simply evaluated using Eq. (54) and one gets the same contribution as in Eq. (p5| 



c 



(1+52)' 



(58) 



Therefore, the general expression for {P^)st depends on the parameter regime considered and it may generally de- 
pend upon the feedback loop high frequency cutoff Wfb and the thermal bath cutoff 137. However, in the common 
experimental situation of a narrow bandwidth around the resonance peak, 7m (1 + 52) < Alj < Wm, and a high Q 
mechanical mode so that the logarithmic correction in Eq. ( p7[ ) can be neglected, the dependence on the frequency 
cutoffs vanishes and one has {P'^)st = {Q'^)st- Therefore, under these conditions, since it is also {QP + PQ)st = 0, 
the stationary state in the presence of the cold damping feedback scheme is an effective thermal state with a mean 
excitation number (n) = 2{Q'^)st — 1/2, where {Q'^)st is given by Eq. (p5|). This effective thermal equilibrium state 
in the presence of cold damping has been already pointed out in [|ll|,^ , within a classical treatment neglecting both 
the back-action and the feedback-induced terms. The present fully quantum analysis shows that cold damping has 
two opposite effects on the effective equilibrium temperature of the mechanical mode: on one hand T is reduced by 
the factor (1 -I- g2)~^, but, on the other hand, the effective temperature is increased by the additional noise terms. 

Let us now consider the optimal conditions for cooling and the cooling limits of the cold damping feedback scheme. 
In the narrow feedback loop bandwidth case, and neglecting the logarithmic correction to (P^)^^^, the stationary 
oscillator energy is given by 



Ust = 2huJrn{Q^ 



4(1+52) 



92 . ^^bT 
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(59) 



This expression coincides with that derived and discussed in pq | , except for the presence of the homodyne detection 
efficiency 77, which was ideally assumed equal to one in p^ ]. The optimal conditions for cooling can be derived in the 
same way as it has been done in [ p^ . The energy Ust is minimized with respect to C, keeping 52 fixed, thereby getting 
Copt — 92I Under these conditions, the stationary oscillator energy becomes 



Vs. 



92 



2 1 



92 



1 



2kBT 1 
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(60) 



showing that, in the ideal limit rj = 1, g2 00 (and therefore C ^ 52 — > 00), also the cold damping scheme is able to 
reach the quantum limit Ust = S.a-'m/2, i.e., it is able to cool the mirror to its quantum ground state, as first pointed 
out in psf . However, differently from the stochastic cooling case of the preceding subsection, the stationary energy 
does not depend on the mechanical quality factor, implying that cooling is easier to achieve using cold damping, 
because the additional condition Q/52 ^ 00 is not necessary in this case. However, cold damping, at variance with 
stochastic cooling feedback, does not yield any nonclassical feature in the steady state. Fig. ^ shows the rescaled 
steady-state energy 2Ust/^i^m versus C plotted for increasing values of 52 (a: 52 — 10, b: 52 — 10"^, c: 52 = 10^, d: 
52 = 10^), with kBT/huJm = 10^ and 77 — 0.8. The figure is essentially indistinguishable from Fig. ^ since, as we have 
seen, the steady states for the two feedback schemes becomes identical for large mechanical quality factors. For high 
gain values, ground state cooling can be achieved also in this case, even with nonunit homodyne detection efficiency. 
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FIG. 5. Rescaled steady-state energy 2Ust/TT,(jJm versus the rescaled input power plotted for different values of g2 (a: 
g-z — 10, b: g-z — 10'', c: g2 = 10^, d; 172 ~ 10^), with kBT/hurn = 10^ and 77 = 0.8. The optimal input power correspondingly 
increases, and for high gain values, ground state cooling can be achieved. 

The ultimate quantum limit of ground state cooling is achieved in both schemes only if both the input power 
and the feedback gain go to infinity. If instead the input power is kept fixed, the effective temperature does not 
monotonically decrease for increasing feedback gain, but, as it can be easily seen from Eqs. ( ^ ) and (p9|), there is 
an optimal feedback gain, giving a minimum steady state energy, generally much greater than the quantum ground 
state energy. The existence of an optimal feedback gain at fixed input power is a consequence of the feedback-induced 
noise term originating from the quantum input noise of the radiation. In a classical treatment neglecting all quantum 
radiation noises, one would have instead erroneously concluded that the oscillator energy can be made arbitrarily 
small, by increasing the feedback gain, and independently of the radiation input power. This is another example of 
the importance of including the radiation quantum noises, showing again that a full quantum treatment is necessary 
to get an exhaustive description of the system dynamics . 

The experimental achievement of ground state cooling via feedback is prohibitive with present day technology. For 
example, the experiments of Refs. ||ll|,|l^ have used feedback gains up to (72 = 40 and an input power corresponding 
to C — 1: a-nd it is certainly difficult to realize in practice the limit of very large gains and input powers. This is 
not surprising, since this would imply the preparation of a mechanical macroscopic degree of freedom in its quantum 
ground state, which is remarkable. The same considerations hold for breaking the standard quantum limit for the 
steady state position fluctuations with the stochastic cooling feedback. 



V. SPECTRAL MEASUREMENTS AND THEIR SENSITIVITY 



Both stochastic cooling and cold damping feedback schemes cool the mirror by overdamping it, thereby strongly 
decreasing its mechanical susceptibility at resonance (see Eqs. (^) and (|5^)). As a consequence, the oscillator does 
not resonantly respond to the thermal noise, yielding in this way an almost complete suppression of the resonance 
peak in the noise power spectrum. Since the effective temperature is proportional to the area below the noise power 
spectrum, this implies cooling. However, the strong reduction of the mechanical susceptibility at resonance means 
that the mirror does not respond not only to the noise but also to any force acting on it. Therefore one expects that 
the SNR of the optomechanical device is not improved by feedback. However, we shall see that this intuitive guess is 
valid only under stationary conditions, and that, at least in the case of an impulsive force, a nonstationary strategy 
can be designed to improve the sensitivity for the detection of a weak classical force. The possibility to use the above 
feedback cooling schemes in a nonstationary way has been first shown in Here we shall reconsider and extend 
the treatment of adopting a general description of nonstationary spectral measurements. 

Spectral measurements are performed whenever the classical force f(t) to detect has a characteristic frequency. 
Since the directly measured quantity is the output homodyne photocurrent Youtit), we define the signal S{oj) as 
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dte-''^' (Youtit)) Fr^it) 



(61) 



where Ft^ (t) is a "filter" function, approximately equal to one in the time interval [0, T„i] in which the spectral 
measurement is performed, and equal to zero otherwise. Using Eq. (pi]), the input-output relation (p^, and the time 
evolution of the position operator Q{t) (Eq. (17a) or (p^)), the signal can be rewritten as 
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where /(w) and Ft^{lu) are the Fourier transforms of the force and of the filter function, respectively, and x(ijj) is 
equal to Xsc{^) or Xcrf(^), according to the feedback scheme considered. 

The noise corresponding to the signal S{uj) will be given by its "variance"; since the signal is zero when /(t) — 0, 
the noise spectrum can be generally written as 
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where the subscript / = means evaluation in the absence of the external force. Using again (|1 1|) , Eqs. (|12|), and the 
input noises correlation functions (|^) and (^3|), the spectral noise can be rewritten as 
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(64) 



where C{t,t') = {Q{t)Q{t') + Q{t')Q{t)) /2 is the symmetrized correlation function of the oscillator position. This 
very general expression of the noise spectrum is nonstationary because it depends upon the nonstationary correlation 
function C{t,t'). The last term in Eq. (BJ) is the shot noise term due to the radiation input noise. 



A. Stationary spectral measurements 

Spectral measurements are usually performed in the stationary case, that is, using a measurement time T,„ much 
larger than the typical oscillator timescalcs. The most significant timescale is the mechanical relaxation time, which 
is in the absence of feedback and [7m(l + 9%)]^^ [i = Ij 2) in the presence of feedback. In the stationary case, the 
oscillator is relaxed to equilibrium and, redefining t' = t + r, the correlation function C{t, t') = C{t,t + T) in Eq. ( |6^ ) 
is replaced by the stationary correlation function Cst{T) = lim4__,oo C{t,t + r). Moreover, for very large T™, one has 
FT,„{t + r) ~ FT,^{t) ~ 1 and, defining the measurement time T™ so that T„i — J dtFT,„{ty, Eq. (Q) assumes the 
form 



N{uj) 



7c 



1/2 



(65) 



where 
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is the stationary position noise spectrum. This noise spectrum can be easily evaluated using the results of the preceding 
section. In fact, using the definition of Cst{T) and the inverse Fourier transform of Eq. (pq), one has 



{Q^)st - lim(Q2(t)) ^c,t(0) 

t — »oo 



(67) 



The position noise spectrum can then be extracted from the stationary mean values derived in the preceding section. 
Using Eq. (|5l|), one has 
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. i . ^ 2'QAU^) + Ti COth 



(68) 



for the cold damping scheme, while the derivation for the stochastic cooling case is less immediate. In fact, using 
Eqs. (|^), (|l]), and (||), one gets 
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Then, using the Fourier transform of Eq. {uM in Eq. fcfl) , one finally gets 
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This position noise spectrum for the stochastic cooling feedback essentially coincides with that already obtained 
in except that in that paper the high temperature limit [cothihio /2kBT) ~ 2kBT /fiuj) is considered and the 
presence of the frequency cutoffs w and w fi, is not taken into account. The noise spectrum in the cold damping case 
of Eq. ( |68| ) instead essentially reproduces the one obtained in [|l8| , with the difference that in Ref . |l8) the homodyne 
detection efficiency 77 is set equal to one, and the feedback and thermal noise cutoff functions have not been explicitely 
considered. The comparison between Eqs. ( |68| ) and (|70| ) shows once again the similarities of the two schemes. The 
only differences lie in the different susceptibilities and in the feedback-induced noise term, which has an additional 
7^/ct;^,j factor in the stochastic cooling case, which is however usually negligible with good mechanical quality factors. 
In fact, it is possible to see that the two noise spectra are practically indistinguishable in a very large parameter 
region. 

The effectively detected position noise spectrum is not given by Eqs. ( |68| ) and (|7^), but one has to add the shot 
noise contribution due to the input noise in the homodyne photocurrent. In fact, using Eq. (|65|), and rescaling it to 
a position spectrum, one has 



Im |Xi(w)|' 



9t 
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■ coth 



Tioj 
2kBT 



1 



4f7C7n 



(71) 



where i = 1 refers to the stochastic cooling case and j = 2 to the cold damping case, The homodyne-detected position 
noise spectrum is actually subject also to cavity filtering, yielding an experimental high frequency cutoff 7c, which 
however does not appear in Eq . (pij) because we have adiabatically eliminated the cavity mode from the beginning. 
Therefore the spectrum of Eq. ([7l|) provides a faithful description of the mirror mode dynamics only for cu < jc', since 
it is usually w,Wfi, > jc, we have not considered the feedback and reservoir cutoff functions in Eq. ([zil), and we 
shall not consider them in the following. The detected noise spectrum has three contributions: the Brownian motion 
term which is independent of the input power p, the shot noise term inversely proportional to p, and the back-action 
term, proportional to p. The main effect of feedback on the spectrum is the modification of the susceptibility due 
to the increase of damping, which is responsible for the suppression and widening of the resonance peak. This peak 
suppression in the noise spectrum has been already predicted and illustrated in and experimentally verified 

for the cold damping case in [^,^. Moreover, the feedback- induced noise term proportional to gf is responsible for 
an increase of the shot noise contribution to the spectrum. For a given feedback gain and frequency, the minimum 
noise is obtained at an intermediate, optimal, power, given by 



■,opt 



and the corresponding value of the minimum displacement noise is 



2iLj„ 



■ coth 



2kBT 



+ + Q-^9l\U^r i^' + kill)- 



(72) 



(73) 



This expression shows that both feedback schemes are able to arbitrarily reduce the displacement noise at resonance. 
In fact, using the fact that Xii'^m) oc in both cases, one has that Nq ^^^{ujm) can be made arbitrarily small 
by increasing the feedback gain. This noise reduction at resonance is similar to that occurring to an oscillator with 
increasing damping, except that in our case, also the feedback-induced noise increases with the gain, and it can be 
kept small only if the input power is correspondingly increased in order to maintain the optimal condition (^2|) . This 
arbitrary reduction of the position noise in a given frequency bandwidth with increasing feedback gain does not hold 
if the input power C is kept fixed. In this latter case, the noise has a frequency-dependent lower bound which cannot 
be overcome by increasing the gain. There is an important difference between the two feedback schemes. In fact, 
it is easy to check from Eq. (^) that in the cold damping case noise reduction takes place only close to resonance, 
and that the noise spectrum is not affected at lower frequencies (for example Nq „j„(ij-' = 0) is not changed by the 



cold damping feedback). In the stochastic cooling case instead, frequency renormalization lo^ 



'Qilm allows to 



reduce position noise even at low frequencies. This reduction of position noise out of resonance, without cold damping 
but with a feedback-induced increase of the mechanical frequency, has been demonstrated experimentally by Cohadon 
et al. in Ref. p|. 
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In the case of stationary spectral measurements also the expression of the signal simplifies. 
Ft,„(uj) ~ 5{uj), and Eq. (p2) assumes the traditional form 



8G/377 



2tt 



In fact, one has 



(74) 



The stationary SNR, TZst{Lo), is now simply obtained dividing the signal of Eq. (|7j) by the noise of Eq. (|65| 

-1/2 
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7mlX«(w)| , 



(75) 



where again i — 1 refers to the stochastic cooling case and i 2 to the cold damping case. It is easy to see that, 
in both cases, feedback always lowers the stationary SNR at any frequency, (except at cj = 0, where the SNR for 
the cold damping case does not depend upon the feedback gain). This is shown in Fig. ^, where the stationary SNR 
in the case of an ideal impulsive force (that is, /(ijj) is a constant) is plotted for three values of the feedback gain. 
The curves refer to both feedback schemes because the two cases i = 1,2 gives always practically indistinguishable 
results, except for very low values of Q. As mentioned at the beginning of the section, this result is not surprising 
because the main effect of feedback is to decrease the mechanical susceptibility at resonance, so that the oscillator is 
less sensitive not only to the noise but also to the signal. Therefore, even though the two feedback schemes are able 
to provide efficient cooling and noise reduction in narrow bandwidths for the mechanical mode, they cannot be used 
to improve the sensitivity of the optomechanical device for stationary measurements. In the next section we shall 
see how cooling via feedback can be used to improve the sensitivity for the detection of impulsive forces, using an 
appropriate nonstationary strategy. 
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FIG. 6. Stationary SNR as a function of frequency in the case of an ideal impulsive force, i.e., f{u!) = const. The full line 
refers to the case with no feedback, the dashed line to the case with 51 = <72 = 10*, and the dotted line to the case with 
gi ~ 92 ~ 10^ (the two feedback schemes give indistinguishable results in these cases). The other parameters are Q = 10^, 
C = 10, ksT/hum = 10^ and r? = 0.8. At a given frequency, the stationary SNR decreases for increasing feedback gain. 



VI. HIGH-SENSITIVE NONSTATIONARY MEASUREMENTS 



The two feedback schemes discussed here achieve noise reduction through a modification of the mechanical sus- 
ceptibility. However, this modification does not translate into a sensitivity improvement because at the same time it 
strongly degrades the detection of the signal. The sensitivity of position measurements would be improved if the oscil- 
lator mode could keep its intrinsic susceptibility, unmodified by feedback, together with the reduced noise achieved by 
the feedback loop. This is obviously impossible in stationary conditions, but a situation very similar to this ideal one 
can be realized in the case of the detection of an impulsive force, that is, with a time duration a much shorter than the 
mechanical relaxation time (in the absence of feedback), a <^ I /7m- In fact, one could use the following nonstationary 
strategy: prepare aX t — Q the mirror mode in the cooled stationary state of Section IV, then suddenly turn off the 
feedback loop and perform the spectral measurement in the presence of the impulsive force for a time Tm, such that 
a <C Tin <SC I/7m- In such a way, the force spectrum is still well reproduced, and the mechanical susceptibility is the 
one without feedback (even though modified by the short measurement time Tm ^ I/7?n)- At the same time, the 
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mechanical mode is far from equilibrium during the whole measurement, and its noise spectrum is different from the 
stationary form of Eq. ([7l|), being mostly determined by the cooled initial state. As long as Tm «C 7m, heating, that 
is, the approach to the hotter equilibrium without feedback, will not affect and increase too much the noise spectrum. 
Therefore, one expects that as long as the measurement time is sufficiently short, the SNR for the detection of the 
impulsive force (which has now to be evaluated using the most general expressions ( |6^ ) and (|6^)) can be significantly 
increased by this nonstationary strategy. 

It is instructive to evaluate explicitely the nonstationary noise spectrum of Eq. ( |6^ ) for the above measurement 
strategy. Let us first consider the cold damping case, which gives more compact expressions. Using Eq. (Ea), one gets 



C{t,t') = K{t)K{t'){Q')st+Xo{t)xo{t'){P')st 



dti 



dt2Xo{ti}xo{h)c{t - t' - ti + t2), 



(76) 



where Xo{t) is the mechanical susceptibility in the absence of feedback (see Eq. ( p4D with = or Eq. ( p3| ) with 
g2 = 0), K(t) is given by Eq. (27) with Xcd replaced by xo, {Q^)st and {P^)st are the stationary values in the presence 
of feedback evaluated in Section IV, and c{t) is the cold damping noise correlation function introduced in Eqs. ( |47| ) 
and (|48|). This nonstationary correlation function has to be inserted in Eq. (|6^). Simple analytical results are obtained 
if we choose the following filter function 



(77) 



{9(t) is the Heavyside step function), satisfying J dtFT^{t)^ — T,n. Using Eq. ( |77| ) and rewriting c{t) in terms of its 
Fourier transform c(w), one gets 



K{uj ~ i/2T„,) {Q^)st + \xo{^ - V2T„0l' {P^)st 



+ \xo{co - l/2Tm)\' 



duj' 



2^ ik 



(78) 



From Eq. (p7[), it is possible to see that K{lu) — (iuj + 7„i)xo('jj)/^m; then, using Eq. ( ^ ) with gd — 0, and the 
high temperature approximation cot]i{huj/2kBT) ~ 2kBT/huj for the Brownian noise, one finally gets the following 
expression for nonstationary noise spectrum for the cold damping feedback 
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{Q^)st 



VTrr 



(79) 



The corresponding noise spectrum for the stochastic cooling case can be obtained in a similar way. Using Eq. (17a), 
one gets 



C{t,t') = KQit)KQ{t'){Q^),t+xo{t)xo{t'){P')st + [xomQit') + KQ{t)xoit')] 



(QP + PQ), 



+ dh dt2Xo{ti)xo{t2)c{t-t' -ti+t2), 



(80) 



where KQ{t) is given by Eq. ( p^ (with Xsc replaced by xo), {Q'^)st, {P'^)st and {QP + PQ)st are the stationary 
values in the presence of stochastic cooling feedback evaluated in Section IV, and we have used the fact that, without 
feedback, ci{t) = and C2{t) = c{t) (see Eqs. ( ^ ) and (^)). Inserting this nonstationary correlation function in 
Eq. (|6^), using Eq. ([77|), the fact that Kq{uli) = {ilu + 7m)xo(^)/'j-'m, and again the high temperature approximation 
for the Brownian noise, one finally gets 
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(81) 



Notice that the two noise spectra ( |79| ) and (81 ) are very similar, the only difference being in the initial stationary values, 
whose explicit expression for the two feedback schemes is given in Section IV. It is also easy to check that the stationary 
noise spectrum corresponding to the situation with no feedback is recovered in the limit of large Tm, as expected, when 
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the terms proportional to jmTm become dominant, and Xo(w — i/2Tm) Xo{^)- In the opposite Hmit of small Tm 
instead, the terms associated to the cooled, initial conditions are important, and since the terms proportional to "fmTm 
are still small, this means having a reduced, nonstationary noise spectrum. This is clearly visible in Fig. 0, where 
the nonstationary noise spectrum, renormalized in order to have a position spectrum, Nq{uj) — N"^ [uj) / ArjC^^mTm: 
is plotted for different values of the measurement time T^, ■ymTm = 10^^ (dotted Hue), jmTm — 10^^ (full line), 
ImTm = 10^"^ (dashed Hue), ^„iT„i = 10~^ (dot-dashed Hue). The resonance peak is significantly suppressed for 
decreasing T^, even if it is simultaneously widened, so that one can even have a slight increase of noise out of 
resonance. This figure is referred to the cold damping feedback scheme, but it is indistinguishable from that obtained 
with the stochastic cooling feedback, using the same parameters (Q = 10*, C — 10, 5i = 32 = 10^, ksT /TiiOm = 10^, 
rj — 0.8). In fact, it can be checked that the two nonstationary noise spectra (79) and ( ^l| ) differ significantly only 
at very low values of the mechanical quality factor (Q < 10^). The effect of the terms depending upon the feedback- 
cooled initial conditions on the nonstationary noise is shown in Fig. 0, where the noise spectrum is plotted for different 
values of the feedback gain at a fixed value of T^. In Fig. ||a, Nq{uj) is plotted at 'jmTm = 10^"^ for 32 = 1 (full 
line), (72 = 10 (dotted line), (72 = 10^ (dashed), (72 = 10'^ (dot-dashed). For this low value of ^mTm, the noise 
terms depending on the initial conditions are dominant, and increasing the feedback gain implies reducing the initial 
variances, and therefore an approximately uniform noise suppression at all frequencies. In Fig. Nq{uj) is instead 
plotted at 'jmTm — 10~^ for 52 = 1 (fuU hnc), 32 = 10 (dotted line), g2 = 10^ (dashed), 52 — 10'^ (dot-dashed). In 
this case, the feedback- gain- independent, stationary terms become important, and the effect of feedback on the noise 
spectrum becomes negligible. Also in this case. Fi g. p| is valid for both stochastic cooling and cold damping schemes. 

It is also possible to check from Eqs. (^9|) and ( ^l[ ) that, similarly to what happens for the stationary case, noise 
does not uniformly decrease for increasing feedback gain if the input power C. is kept fixed, but there is an optimal 
feedback gain, minimizing the noise at a given frequency and input power. 
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FIG. 7. Nonstationary noise spectrum 



A'' (ajm)/477("7,„r,„ for different values of the measurement time, 
TmTm = 10" (dotted line), 7^Tm = 10"^ (full line), 7mTm = 10"^ (dashed line), 7™r„ = 10"* (dot-dashed line). The 
figure refers to the cold damping feedback scheme, but the curves are indistinguishable from that obtained with the stochastic 
cooling feedback, using the same parameters, Q = 10*, — 10, Qi ~ 92 = 10^, kBT/hujm. ~ 10^, ry — 0.8. 
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FIG. 8. Nonstationary noise spectrum Nq{lo) for different values of the feedback gain, g2 ~ 1 (full line), g2 = 10 (dotted 
line), (?2 = 10^ (dashed), (?2 = 10^ (dot-dashed), with fixed measurement time, "fmTm = 10~^ (a), and 'jmTm ~ 10"^ (b). (a) 
corresponds to a strongly nonstationary condition, in which the noise is significantly suppressed, thanks to the cooled initial 
condition. In (b) the stationary terms becomes important and the noise reduction due to feedback cooling is less significant. The 
figure refers to the cold damping feedback scheme, but the curves are indistinguishable from that obtained with the stochastic 
cooling feedback, using the same parameters, Q = 10*, = 10, kBT/hujm ~ 10^, rj — 0.8. 



The significant noise reduction attainable at short measurement times jmTm ^ 1 is not only due to the feedback- 
cooled initial conditions, but it is also caused by the effective reduction of the mechanical susceptibility given by the 
short measurement time, Xq{i-u) — > Xo(^ ~ V^Tm)- This lowered susceptibility yields a simultaneous reduction of the 
signal at small measurement times 7m Tm ^ 1, and therefore the behavior of the nonstationary SNR may be nontrivial. 
However, one expects that impulsive forces at least can be satisfactorily detected using a short measurement time, 
because the noise can be kept very small and the corresponding sensitivity increased. Let us check this fact considering 
the case of the impulsive force 

fit) = /o exp [~it ~ ti) VSa^] cos (uft) , (82) 

where a is the force duration, ti its "arrival time", and ujf its carrier frequency. The corresponding SNR is obtained 
dividing the signal of Eq. (|6l]), evaluated with Eq. (77), by the nonstationary noise spectra of Eqs. ( [79[ ) and (^ij), and it 



is shown in Figs. || and |T^. As anticipated, the sensitivity of the optomechanical device is improved using feedback in a 
nonstationary way. In Fig. ^, the spectral SNR, TZ{uj), is plotted for different values of feedback gain and measurement 
time (as in the previous curves, the figures well describe both feedback schemes, because they give indistinguishable 
results for TZ(lo) in the physically relevant parameter region). The full line refers to = <72 = .9 = 2 • 10"^ and 
ImTm = 10"'^, the dashed line to the situation with no feedback and the same measurement time, g = and 
ImTm = 10~^; finally the dotted line refers to a "standard" measurement, that is, no feedback and a stationary 
measurement, with a long measurement time, jmTm = 10. The proposed nonstationary measurement scheme, "cool 
and measure", gives the highest sensitivity. This is confirmed also by Fig. |l^, where the SNR at resonance, TZ{ujm), 
when feedback cooling is used with g = 2-10'' (full line), and without feedback cooling (dotted line), is plotted as 
a function of the rescaled measurement time "fmTm- The preparation of the mirror in the cooled initial state yields 
a better sensitivity for any measurement time. As expected, the SNR in the presence of feedback approaches that 
without feedback in the stationary limit 7m Tm 3> 1, when the effect of the initial cooling becomes irrelevant. Both 
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Fig. ^ and |l^ refer to a resonant {uf = Um) impulsive force with 7^(7 = 10 ^ and j„iti = 3-10 ^, while the other 
parameters are Q — 10^, C = 10, 77 = 0.8, kBT/hujm — 10^. 
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FIG. 9. Spectrum of the nonstationary SNR, TZ{u>), with and without feedback cooling of the initial state. The full line 
refers to a nonstationary measurement, 7mTm = 10~^, in the presence of feedback, g — 2 ■ 10"^ (the two feedback schemes 
give indistinguishable curves); the dashed line refers to the no-feedback case, and with the same, short, measurement time 
7m7m = 10^'^. Finally, the dotted line refers to a "standard measurement", without feedback, and in the stationary limit 
7mT™ = 10. The other parameters are cjf = oj™, 7™cr = 10"*, 7„ti = 3 • 10"*, Q = 10^ C = 10, r; = 0.8, kBT/hujm = 10^. 
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FIG. 10. Nonstationary SNR at resonance, TZ{uJrn), with and without feedback cooling of the initial state, plotted as a 
function of the rescaled measurement time 7mTm. The full line refers to the case with feedback-cooled initial conditions 
{g = 2 ■ 10^, the two feedback schemes give indistinguishable curves). The dotted line refers to the no- feedback case, g = 0. 
The other parameters are the same as in Fig. ^ 

The proposed nonstationary strategy can be straightforwardly applied whenever the "arrival time" ti of the impul- 
sive force is known: feedback has to be turned off just before the arrival of the force. However, the scheme can be easily 
adapted also to the case of an impulsive force with an unknown arrival time, as for example, that of a gravitational 
wave passing through an interferometer. In this case it is convenient to repeat the process many times, i.e., subject 
the oscillator to cooling-heating cycles. Feedback is turned off for a time T,„ during which the spectral measurement is 
performed and the oscillator starts heating up. Then feedback is turned on and the oscillator is cooled, and then the 
process is iterated. This cyclic cooling strategy improves the sensitivity of gravitational wave detection provided that 
the cooling time Tcooi, which is of the order of 1/ [7m (1 + gi)], is much smaller than T,„, which is verified at sufficiently 
large gains. Cyclic cooling has been proposed, in a qualitative way, to cool the violin modes of a gravitational waves 
interferometer in , and its capability of improving the high-sensitive detection of impulsive forces has been first 
shown in [T^ ]. In the case of a random, uniformly distributed, arrival time ti and in the impulsive limit a <C Tm, the 
performance of the cyclic cooling scheme is well characterized by a time averaged SNR, i.e., 
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(83) 
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where TZ{uj,ti) is the nonstationary SNR at a given force arrival time ti discussed in this section, and TZ{uj,ti)cooi is 
the nonstationary SNR one has during the cooling cycle, which means with feedback turned on and with uncooled 
initial conditions. It is easy to understand that TZ{u!,ti)cooi ^ TZ{uj,ti), and, since it is also Tcooi ^ ?m, the second 
term in Eq. |3|) can be neglected, so that [ p^ , 

-i-rn ^ -i- cool J 

This time-averaged SNR can be significantly improved by cyclic cooling, as it is shown in Fig. |ll|, where {TZ{uj)) is 
plotted both with and without feedback. The full line describes the time-averaged SNR subject to cyclic feedback- 
cooling with g = 2 ■ 10^, ^mTm = lO"'^, and Tcooi = IQ-'^Tm. In the absence of feedback, in the case of an impulsive 
force with unknown arrival time and duration a, the best strategy is to perform repeated measurements of duration 
Tm without any cooling stage. The measurement time Tm can be optimized considering that it has to be longer than 
(T, and at the same time it has not to be too long, in order to have a good SNR (see the dotted line in Fig. |l^). In 
this case, the time-averaged SNR can be written as 

(7^o(t^)) -7^ [ dt{llo{uM). (85) 

where TZo{u!,ti) is the SNR evaluated for g = 0. The dashed line in Fig. |ll| refers to this case without feedback, and 
with 'jmTm — 10"^. The other parameter values are the same as in Figs. ^ and |l^ and in this case, cyclic cooling 
provides an improvement at resonance by a factor 16 with respect to the case with no feedback. As suggested in 
Ref. [|l5|, one could use nonstationary cyclic feedback to cool the violin modes in gravitational- wave interferometers, 
which have sharp resonances within the detection band. One expects that single gravitational bursts, having a duration 
smaller than the cooling cycle period, could be detected in this way. 
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FIG. 11. Time averaged spectral SNR with and without cyclic cooling. The full line refers to cyclic cooling with ^mTm = 10" , 
g — 2 ■ 10^, and Tcooi = W~^Tm (the two feedback schemes give indistinguishable curves). The dashed line refers to the 
no- feedback case, with the same measurement time "fmTm = 10"'' (see Eq. (|8^)). The other parameters are LUf = uj,-a, 
7™cr = 10"*, Q = 10^ ( = 10, ri = 0.8, kBT/hiUm = 10'\ 



VII. CONCLUSIONS 



We have studied how quantum feedback schemes can be used to reduce thermal noise and improve the sensitivity 
of optomechanical devices. We have analysed in detail the stochastic cooling scheme introduced in Ref. and the 
cold damping scheme experimentally implemented in Ref. We have seen that the two schemes are physically 

analogous, even though they show some differences. In both cases, the main effect of feedback is the increase of 
mechanical damping, accompanied by the introduction of a controllable, measurement-induced, noise. The increase of 
damping means reduction of the susceptibility at resonance, and the consequent suppression of the resonance peak in 
the noise spectrum. Stochastic cooling feedback differs form cold damping in the fact that it has the supplementary 
effect of increasing the mechanical frequency. This means that, while cold damping achieves thermal noise reduction 
only around resonance, stochastic cooling is able to reduce noise even at very low frequencies, out of resonance. We 
have also shown that both schemes are able to achieve the ultimate quantum limit of ground state cooling (see also 
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|l8[ for the cold damping case). For both feedback schemes, ground state cooling is reached in the limit of very large 
feedback gain, ideal homodyne detection, and very large input power. In the stochastic cooling case, however, also 
the additional condition of very large mechanical quality factor is needed (see also |Q ) , so that cooling is much more 
easily achieved in the cold damping case. In the limit of very large gain and input power, but with fixed mechanical 
quality factor, stochastic cooling feedback is instead able to achieve steady state position squeezing, that is, one can 
beat the standard quantum Hmit {Qlf) < 1/4. Finally stochastic cooling is also able to produce stationary contractive 
states ]36[ |. Reaching these quantum limits in optomechanical sytems is experimentally very difficult but it would be 
extremely important, because it would be a genuine manifestation of quantum mechanics for a macroscopic mechanical 
degree of freedom. 

We have also analysed the sensitivity of the optomechanical device in the case of position spectral measurements for 
the detection of weak forces. Even though both feedback schemes are not able to improve the sensitivity of stationary 
measurements, we have shown how feedback can be used in a nonstationary way in order to increase of the SNR 
in the case of impulsive forces. If the arrival time of the classical force is known, one has to keep the mirror mode 
cooled by feedback, and then turn off the feedback just before the arrival of the force. The mirror therefore responds 
to the force with its intrinsic susceptibility, not suppressed by the feedback, and with a nonstationary noise, reduced 
by the feedback. The SNR is increased as long as the measurement time Tm is longer than the force duration tr, 
but much smaller than the mechanical relaxation time, that is, a Tm ^ l/7m- This nonstationary strategy can 
be well adapted to the case of a force with an unknown arrival time, as for example, gravitational waves. In this 
case, the cooling and measurement steps has to be cyclically repeated, and the performance of cyclic cooling can be 
characterized by a SNR averaged over the force arrival time. This time-averaged SNR can be significantly improved 
by cyclic cooling, thanks also to the fact that the cooling time can be made very small using very large feedback 
gains 5, because it is Tcooi — [7m(l + g)] ^ ■ Differently from ground state cooling, the experimental implementation 
of these nonstationary strategy is feasible with current technology, and it may be useful not only for optomechanical 
devices, but also for microelectromechanical systems. 
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